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(57) ABSTRACT 

The present invention relates to a space-invariant independent 
component analysis and electronic nose for detection of 
selective chemicals in an unknown environment, and more 
specifically, an approach to analysis of sensor responses to 
mixtures of unknown chemicals by an electronic nose in an 
open and changing environment. It is intended to fill the gap 
between an alarm, which has little or no ability to distinguish 
among chemical compounds causing a response, and an ana- 
lytical instrument, which can distinguish all compounds 
present but with no real-time or continuous event monitoring 
ability. 
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Correlation Between The original And Separated Sources of Chemicals 
3 & 10 (Bold entries indicate source chemical with highest correlation) 
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Correlation Between The Original And Separated Sources of Chemicals 
2 & 5 (Bold entries indicate source chemical with highest correlation) 
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Projection Mean & Standard Deviation for Chemicals 1 & 7 for each Source 
Clearly Identify Proper Source (Bold entries indicate highest correlation) 
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METHOD FOR DETECTION OF SELECTED 
CHEMICALS IN AN OPEN ENVIRONMENT 

PRIORITY CLAIM 

5 

The present application is a non-provisional utility appli- 
cation, claiming benefit of priority of U.S. Provisional Patent 
Application No. 60/715,353, filed Sep. 7, 2005, titled “Smart 
Enose for Chemical Detection in the Open Environment.” 

to 

STATEMENT OF GOVERNMENT INTEREST 

The invention described herein was made in the perfor- 
mance of work under a NASA contract and is subject to the 
provisions of Public Law 96-517 (35 USC 202) in which the 15 
Contractor has elected to retain title. 

BACKGROUND OF THE INVENTION 

(1) Technical Field 20 

The present invention relates to a space-invariant indepen- 
dent component analysis and electronic nose for detection of 
selective chemicals in an unknown environment, and more 
specifically to an approach to analysis of sensor responses to 
mixtures of unknown chemicals by an electronic nose in an 25 
open and changing environment. 

(2) Background 

The need for low-power, miniature sensor devices that can 
monitor air quality in an enclosed space with multi-com- 
pound capability and minimum human operation led to the 30 
development of a polymer-carbon composite based electronic 
nose (ENose) at NASA’s Jet Propulsion Laboratory (JPL). 
The sensor array in the JPL ENose consists of 32 conducto- 
metric sensors made from insulating polymer films loaded 
with carbon. In its current design, it has the capability to 35 
detect 1 0 co mm on contaminants which may be released into 
the recirculated breathing air of a space shuttle or space 
station released from a spill or a leak; target concentrations 
are based on the 1-hour Spacecraft Maximum Allowable 
Concentrations (SMAC) set by NASA, depicted in FIG. 1, 40 
and are in the parts-per-million (ppm) range. The ENose was 
intended to fill the gap between an alarm, which has little or 
no ability to distinguish among chemical compounds causing 
a response, and an analytical instrument, which can distin- 
guish all compounds present but with no real-time or continu- 45 
ous event monitoring ability. 

As in other array -based sensor devices, the individual sen- 
sor films of the ENose are not specific to any one analyte; it is 
in the use of an array of different sensor films that gases or gas 
mixtures can be uniquely identified by the pattern of mea- 50 
sured response. The response pattern requires software analy- 
sis to deconvolute gas compounds and their concentrations. 

An example sensor set is shown on the left in FIG. 2, with 
the complete assembled device shown on the right. 

What is needed is a method of detection for selective 55 
chemicals as a result of leaks or spills of specific compounds. 

It has been shown in analysis of samples taken from space 
shuttle flights that, in general, air is kept clean by the air 
revitalization system and contaminants are present at levels 
significantly lower than the SMACs; the ENose has been 60 
developed to detect target compounds released suddenly into 
the breathing environment. A leak or a spill of a solvent or 
other compound would be an unusual event. 

What is needed is an approach to analysis of sensor 
responses to mixtures of chemical compounds so that use of 65 
the ENose may by extended to detect chemical compounds in 
an open and changing environment, such as a building or a 


2 

geographical area where air exchange is not controlled and 
limited. In an open environment, the collected sensory data 
will be comprised of a mixing between unknown chemicals 
with unknown mixing levels (coefficient) between them. The 
identification of chemical compounds among these mixing 
chemicals is a challenge for real world applications. 

To determine whether a chemical compound exists in the 
an environment, one of the most well-known techniques is to 
recover the original chemicals. When done, the detection can 
be an easy step by determining the minimum phase between 
the predicted original reactants and the target chemicals. A 
more sophisticated method is to use a neural network 
approach, which canbe employed to capture the target chemi- 
cals in various conditions, e.g., concentration levels through 
the parameterized weight set. Then, the strongest correlation 
between parameterized weight and the predicted original can 
be used to identify the intended chemical. 

Recently, Independent Component Analysis (ICA) has 
proven effective to not only de-correlate second order statis- 
tics of the signals but also reduce higher order statistical 
dependencies. ICA transforms an observed signal vector into 
a set of signals that are as statistically independent as pos- 
sible. Theoretically, ICA is an information-theoretic 
approach, which exploits concepts from information theory 
such as entropy and mutual information. 

The ICA roots in the early work of Herault and Jutten who 
first introduced an adaptive algorithm in a simple feedback 
architecture that was able to separate several unknown inde- 
pendent sources. ICA was further developed, and recent 
improvements used natural gradient descent based on the 
Riemannian metric tensor to optimize the curvature of a par- 
ticular manifold in n dimensional space. This technique is 
employed to apply to the Infomax to simplify the learning rule 
used here. ICA has applications for feature extraction in 
speech recognition systems, in communication systems, in 
medical signal processing, and in image processing. 

Therefore, what is needed is an ICA method for detection 
of selective chemicals in an unknown environment using an 
electronic nose. 

SUMMARY OF THE INVENTION 

The present invention relates to a space-invariant indepen- 
dent component analysis and electronic nose for detection of 
selective chemicals in an unknown environment, and more 
specifically, to an approach to analysis of sensor responses to 
mixtures of unknown chemicals by an electronic nose in an 
open and changing environment. 

In one aspect of the present invention, a method for detec- 
tion of chemicals in an open environment comprises the acts 
of obtaining sampled data by sampling data from a database 
containing air samples of an unknown open environment; 
utilizing an electronic nose sensor array containing a plurality 
of sensors evaluating the sampled air over a plurality of time 
intervals; partitioning the sampled data with respect to each 
time interval into a subset; sampling a subset of the sampled 
data with respect to time for each sensor; calculating operat- 
ing points of the subset of the sampled data as the averaged 
value (e.g., mean value) of the subset of the sampled data; 
linearizing the subset of the sampled data using the operating 
points; determining a data distribution corresponding to the 
linearized subset of the sampled data; performing an indepen- 
dent component analysis (ICA) on the linearized subset of the 
sampled data by using the data distribution, wherein the ICA 
generates independent component vectors (recovered 
sources) representing the subset of the sampled data; sam- 
pling known chemical elements representing spectral vectors 
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from a database; iteratively performing an act of projecting 
the independent component vectors (recovered source) rep- 
resentative of the subset of the sampled data onto the inde- 
pendent component vectors of the known chemical elements, 
wherein the step of projecting further comprises the steps of: 5 
finding a correlation between the independent component 
vector of a subset of the sampled data and the spectral vectors 
of the known chemical elements; finding a maximum corre- 
lation point on the correlation and a known chemical element 
corresponding to the maximum correlation point; creating a to 
list of chemical elements detected on the subset of the 
sampled data, wherein a known chemical element corre- 
sponding to the maximum correlation point is added to the list 
of chemical elements detected; outputting the list of chemical 
elements detected on the subset of the sampled data. 15 

BRIEF DESCRIPTION OF THE DRAWINGS 

The objects, features and advantages of the present inven- 
tion will be apparent from the following detailed descriptions 20 
of the disclosed aspects of the invention in conjunction with 
reference to the following drawings, where: 

FIG. lisa table depicting the 1 -hour Spacecraft Maximum 
Allowable Concentrations (SMAC) of chemicals as deter- 
mined by NASA; 25 

FIG. 2 is a photograph of a sensor set of an electronic nose 
(Enose) and assembled Enose device; 

FIG. 3 is a diagram of the system model and setup topology 
of an Enose array; 

FIG. 4 is a diagram depicting the difference between a 30 
time-invariant independent component analysis (ICA) of the 
prior art and the space-invariant ICA architecture of the 
present invention; 

FIG. 5 is a graphical illustration of the spectra of eleven 
chemicals used for comparison with mixtures of chemicals 35 
from an unknown source; 

FIG. 6 is a table depicting the four sets of data used to detect 
selective chemicals from chemical mixtures; 

FIG. 7A is a graphical representation of twelve mixtures of 
chemicals 1 and 7 used in Experiment 1 ; 40 

FIG. 7B is a graphical representation of separated signal 
sources of chemicals 1 and 7 by space-invariant ICA of 
twelve inputs; 

FIG. 7C is a graphical representation of the average of the 
single chemical source 1 and 7; 45 

FIG. 8 is a table depicting a projection of the separated 
chemical sources 1 and 7 on the original chemical; 

FIG. 9A is a graphical representation of 12 mixtures of 
chemicals 3 and 10 used in Experiment 1; 

FIG. 9B is a graphical representation of separated signal 50 
sources of chemicals 3 and 10 by space-invariant ICA of 12 
inputs; 

FIG. 9C is a graphical representation of the average of the 
single chemical source 3 and 10; 

FIG. 10 is a table depicting a projection of the separated 55 
chemical sources 3 and 10 on the original chemical; 

FIG. 11 A is a graphical representation of 12 mixtures of 
chemicals 2 and 13 used in Experiment 1 ; 

FIG. 1 IB is a graphical representation of separated signal 
sources of chemicals 2 and 13 by space-invariant ICA of 12 60 
inputs; 

FIG. 11C is a graphical representation of the average of the 
single chemical source 2 and 13; 

FIG. 12 is a table depicting a projection of the separated 
chemical sources 2 and 5 on the original chemical; 65 

FIG. 13 A is a graphical representation of 12 mixtures of 
chemicals 2 and 5 used in Experiment 1 ; 
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FIG. 13B is a graphical representation of separated signal 
sources of chemicals 2 and 5 by space-invariant ICA of 12 
inputs; 

FIG. 13C is a graphical representation of the average of the 
single chemical source 2 and 5; 

FIG. 14 is a table depicting a projection of the separated 
chemical sources 2 and 5 on the original chemical; and 

FIG. 15 is a table depicting the mean and standard devia- 
tion of its projection of separated sources of chemical 1 and 
chemical 7 on each single chemical source, demonstrating the 
successful identification of the separated chemical sources. 

DETAILED DESCRIPTION 

The present invention relates to a space-invariant indepen- 
dent component analysis and electronic nose for detection of 
selective chemicals in an unknown environment, and more 
specifically, to an approach to analysis of sensor responses to 
mixtures of unknown chemicals by an electronic nose in an 
open and changing environment. 

The following description, taken in conjunction with the 
referenced drawings, is presented to enable one of ordinary 
skill in the art to make and use the invention and to incorpo- 
rate it in the context of particular applications. Various modi- 
fications, as well as a variety of uses in different applications, 
will be readily apparent to those skilled in the art, and the 
general principles, defined herein, may be applied to a wide 
range of embodiments. Thus, the present invention is not 
intended to be limited to the embodiments presented, but is to 
be accorded the widest scope consistent with the principles 
and novel features disclosed herein. Furthermore, it should be 
noted that unless explicitly stated otherwise, the figures 
included herein are illustrated diagrammatically and without 
any specific scale, as they are provided as qualitative illustra- 
tions of the concept of the present invention. 

In the following detailed description, numerous specific 
details are set forth in order to provide a more thorough 
understanding of the present invention. However, it will be 
apparent to one skilled in the art that the present invention 
may be practiced without necessarily being limited to these 
specific details. In other instances, well-known structures and 
devices are shown in block diagram form, rather than in 
detail, in order to avoid obscuring the present invention. 

The reader’s attention is directed to all papers and docu- 
ments that are filed concurrently with this specification and 
are open to public inspection with this specification, and the 
contents of all such papers and documents are incorporated 
herein by reference. All the features disclosed in this specifi- 
cation, (including any accompanying claims, abstract, and 
drawings) may be replaced by alternative features serving the 
same, equivalent or similar purpose, unless expressly stated 
otherwise. Thus, unless expressly stated otherwise, each fea- 
ture disclosed is one example only of a generic series of 
equivalent or similar features. 

Furthermore, any element in a claim that does not explic- 
itly state “means for” performing a specified function, or 
“step for” performing a specific function, is not to be inter- 
preted as a “means” or “step” clause as specified in 35 U.S.C. 
Section 112, Paragraph 6. In particular, the use of “step of’ or 
“act of’ in the claims herein is not intended to invoke the 
provisions of 35 U.S.C. 112, Paragraph 6. 

The present invention provides a space-invariant topology 
to enable an Independent Component Analysis (ICA) to solve 
chemical detection from two unknown mixing chemical 
sources. The specific system architecture enables the space- 
invariant ICA to achieve a robust independent signal source 
separation in an unknown open environment. Additionally, 
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the system architecture provides an optimal topology for 
hardware-friendly implementation and an effective architec- 
ture for reduced data source requirements. 

The space-invariant topology along with the techniques of 
maximum entropy information and natural gradient descent 5 
demonstrates that it is effective to separate the two mixed 
unknown chemical sources with unknown mixing levels to 
the array of two original sources under insufficient sampled 
data. From separated sources, they can be identified by pro- 
jecting them on 1 1 known chemical sources to find the best 10 
match for detection. 

Simulations have shown that 1 00% correct detection could 
be achieved under two cases: a) under-completed case where 
the number of input mixtures is larger than number of original 
chemical sources; and b) regular case where the number of 15 
input mixtures is the same as the number of original chemical 
sources. The time-invariant topology approach may face 
obstacles with an over-complete case, insufficient data and 
cumbersome topology. 

20 

Four sets of unknown paired mixture sources are collected 
via an electronic nose, which, in one non-limiting example, is 
a JPL 16-ENose sensor array (made at the Jet Propulsion 
Laboratory, Pasadena, Calif.). The Enose collects the mixture 
sources in the unknown environment with, at most, 12 

25 

samples’ data collected. Per time-invariant aspect, this 
appears to be an over-complete case in ICA where the number 
of outputs (32) is larger than the number of inputs (16). 

Technical Approach 

In one non-limiting example, the approximation of the set 30 
up topology of the electronic nose array is shown in FIG. 1. 
The collected sensing data x z (t) consists of changes in elec- 
trical resistance corresponding to sensor response to the 
unknown mixture of chemical sources s. and their densities 
(or concentrations) at the time t. 35 

Due to the small separation between the sensors them- 
selves, the input of each sensor resistance is assumed to be 
uniquely distributed. The sensory data can be modeled as 
follows: 


*«(0 = Yj a J s jW 


/= l -k 


( 1 ) 


45 


Where f. is the unknown non-linear activation, c^. is the 
unknown mixing coefficient of chemical source j, and i is the 
index of sensor number and N as number of sensor and s^. is 
original source. 50 

From equation (1), using the first order of a Taylor expan- 
sion, it is rewritten: 


*«■ = a > + Yj &j s j 

j = i 


N \ N 


Where a, = / I y atf + £ 

\j = 1 J J = 1 ' 


ajsTj 


df 

and Ay = aj 

dsj & 


(2) 55 


60 


(3) 


And s °j is an operating point of the source s,-. For each 65 
sampling data point in time for the same sensor i, x z - fluctuates 
around its operating point a z and it can be considered as a 


common bias for x z (t) with te[t, t+kAt]. From this argument, 
Equation (2) can be simplified to: 


x { -a { 


Y = 


x 2 -a 2 


x k a k 


Y = rs 

Where Y is unbiased mixture data, T is the unknown mix- 
ing matrix and S is the chemical source signal. 

The learning rule based on the maximum entropy algo- 
rithm is given by: 

y=g(u)=g(Wx); (5) 

where g is a non linear function e.g the logistic function or 
hyperbolic tangent function. The update weight can be cal- 
culated as: 

A W=Vr T +<&{u)x?; (6) 

where W _r is an inverse transport of the NXN weight 
matrix W, x r is a mixing input vector (observed vector), and 


$(«) = [<Pi («/)] = 


' dy’i 


(7) 


To simplify equation (5) using a natural gradient descent 
technique, the learning rule can be: 

A W=(I+®(u)u T )W (8) 

With u=Wx 

Space-Invariant Architecture 

The most common ICA approach is that the number of 
variables and the number of sources are the same. However, in 
this study there are two obstacles: 1) there are 12 or less 
samples (mixing chemical compounds) from each sensor and 
the total number of sensor is 1 6 and they do not have sufficient 
data set (at least 16 data samples required); and 2) the number 
of variables is 16 as number of sensors while the number of 
compounds in a mixture is 2 and it is considered over com- 
plete case. 

For the time-invariant approach, the data that will flow 
orthogonal 402 with the time-invariant direction as shown in 
FIG. 4 and will require 32 outputs (16 channels for each 
chemical compound). The topology is 16 inputs, 32 outputs 
and 12 or less sample data which may not be a solvable 
problem. 

Instead of using a time-invariant approach, the space-in- 
variant approach 404 allows for more data points and enables 
the square mapping matrix (the dimension of mixing sources 
and sensors are the same). This approach is feasible due to the 
mathematical model based on equations (2)- (4). The archi- 
tectures are shown in FIG. 4. In FIG. 4, the unbiased input 
Y z (t) (i=l-k) is based on temporal mixture data and the sen- 
sory data are spatially invariant. 

From laboratory set up, a set of single spectra of 1 1 chemi- 
cals using 1 6 elements is collected in the ENose sensor array; 
it is averaged and shown in FIG. 5. A database may also 
contain the sampled data of air samples from an unknown 
environment to be used in the calculations. 

In this embodiment, the number of sensors used is 16 and 
mixing chemical sources is 2; four sets of data will be exam- 
ined, as shown in table of FIG. 6. 
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Simulation Results 

Based on the data available provided in FIG. 6, two experi- 
ments were performed: 

Experiment 1 : Over-complete case: in this case, a complete 
sample data (see column 4 in FIG. 6) is used as input to 
the network and the output size is 2, as two original 
sources recovered. 

Experiment 2: Squared case: this is a straight forward case 
with 2 inputs and 2 outputs when the data was rearranged 
so that each input is from the same sensor with non- 
overlapped consecutive sampling times (t+i*At and 
t+(i+l)*At). 

Experiment 1 

In this case, the four subsets of the sampled data in FIG. 6 
are studied and all data vectors available are used (the maxi- 
mum number of data vectors is 12), which is less than the 
number of sensors (16). 

Data Set 1 

For data set 1, there are 12 mixtures of chemicals 1 and 7 
and the data are shown in FIG. 7A. 

FIG. 7B shows the separated signal sources 1 702 and 7 704 
via space-invariant ICA of 12 inputs. FIG. 7C shows the 
original chemical sources 1 706 and 7 708 by averaging 
techniques. 

Using the space-invariant ICA approach, the recovered 
signal sources (chemical 1 and 7 sources) are shown in FIG. 
7B and the average of the single chemical source 1 and 7 are 
shown in FIG. 1C. 

The spectrum of signal 702 in FIG. 7B is most closely 
matched to the signal 706 in FIG. 1C. Similarly the spectrum 
of signal 704 in FIG. 7B is most closely matched to the signal 
708 in FIG. 1C. 

To confirm its performance, the separated sources 1 and 7 
were projected by ICA technique on the known 1 1 chemical 
sources shown in FIG. 5, the results are provided in the table 
in FIG. 8. 

As can be inferred from FIG. 8, the single source of chemi- 
cal 1 has the greatest overlap with the separated source, 
labeled separated chemical 1. Similarly, single source chemi- 
cal 7 has the greatest overlap with separated chemical 7 from 
the mixture shown in FIG. 7A. 

Data Set 2 

Data set 2, a mixture of chemicals 3 and 10, is plotted in 
FIG. 9 A below. 

For the mixture from data set 2, the performance of space- 
invariant ICA has demonstrated its effective capability to 
separate the mixture of chemicals 3 and 10 as shown in FIG. 
9B as compared with the original sources in FIG. 9C. 

Table IV in FIG. 10 clearly indicates that the maximum 
sensor value in the column for separated chemical 3 was for 
chemical 3 while the maximum sensor value in the separated 
chemical 10 column was for chemical 10. 

Data Set 3 

Data set 3, a mixture of chemicals 2 and 13, is plotted in 
FIG. 11A. 

For the mixture of chemicals 2 and 13, space-invariant ICA 
has demonstrated its capability to separate the mixture of 
chemical 2 and 13 as shown in FIG. 11B as compared with the 
original sources in FIG. 11C. 

FIG. 12 also clearly indicates that the separated chemical 
with the highest sensor value was chemical 2 among the 
separated chemical 2 data while chemical 13 captured the 
highest value among the data for separated chemical 13. 


Data Set 4 

Data set 4, a mixture of chemicals 2 and 5, is plotted in FIG. 
13A. 

For mixture from data set 4 shown in FIG. 13A, space- 
5 invariant ICA has confirmed its powerful capability to sepa- 
rate the mixture of chemicals 2 and 5 as shown in FIG. 13B as 
compared with the original sources in FIG. 13 C. 

Table VI of FIG. 14 shows the correlation between the 
original and separated source of chemical 2 and 5. 
to Experiment 2 

In this experiment, data set 1 (1 6x12) is paired in columns 
to obtain the data set 96x2. From this conversion, data values 
in a single row are the data from the same sensor with con- 
secutive sampling times (t+i*At and t+(i+l)*At); this new 
15 data set allows the same number of mixing sources and of 
original sources. 

Using this new data set, space-invariant ICA has produced 
the results that were validated with the 1 1 known chemicals. 
Classification was 100% correct, based on the projection on 
20 11 classes shown in FIG. 5. 

To simplify the results, the mean and standard deviation of 
its projection were tabulated (separated sources of chemical 1 
and chemical 7) on each single chemical source, and the 
results are summarized in the table in FIG. 15. 

25 FIG. 15 is compact information to show that the separation 
sources from the mixture have successfully identified the 
original chemical sources. 

Appendix A, incorporated herein by reference, further 
includes a series of flow chart diagrams providing the order in 
30 which data is processed and analyzed in the space-invariant 
independent component analysis. 

Discussion 

To separate two (2) mixing sources from a sixteen (16) 
element sensory data array, known as the over-complete case, 
35 poses a challenges for mathematical model and network 
topology. The non-overlapped paire-wise (i.e. sensor i and 
sensor i+1) or overlapped pair-wise (i.e. (sensor i and sensor 
i+1) and (sensor i+1 and sensor i+2) so on) may face cum- 
bersome and ineffective techniques. 

40 As shown above, the mathematical model has demon- 
strated space-invariant ICA to be an effective topology to 
overcome insufficient data samples and the over-complete 
case. Moreover, the chemical data itself is fuzzy and incon- 
sistent. The simulation demonstrated that the chemical source 
45 separation problem can be solved effectively with complete 
time sampling data (k=12) (under complete case) and two 
consecutive sampling data (k=2). Optimal topology may 
require a model of noise in order to determine the size of the 
sampling input. Moreover, space-invariant ICA governed by 
50 equation (4) is only valid when the sampling time is suffi- 
ciently small. Hence, the sampling time also plays an impor- 
tant role to ensure that the model approach holds. 

Conclusions 

A successful mathematical model to enable the space- 
55 invariant ICA topology from which Infomax and natural gra- 
dient descent technique can be applied has been demon- 
strated, and simulation has confirmed that the modeling is 
effective and sufficient to perform chemical source separation 
to enable the smart ENose to detect mixtures of chemicals in 
60 an open environment. 

What is claimed is: 

1. A computer implemented method for detection of 
chemicals in an open environment, the method comprising an 
65 act of causing a computer to perform operations of: 

obtaining sampled data from a database containing data of 
air samples of an unknown open environment, the data- 
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base created using a sensor array containing a plurality 
of sensors to obtain the data of air samples over a plu- 
rality of time intervals; 

sampling a subset of the sampled data with respect to time 
for each sensor, wherein each subset includes a plurality 5 
of elements; 

calculating operating points of each subset of the sampled 
data as an averaged value of each subset of the sampled 
data; 

linearizing each subset of the sampled data using the oper- 1 c 
ating points by subtracting the averaged value from each 
element in each corresponding subset of the sampled 
data; 

determining a data distribution corresponding to each lin- 
earized subset of the sampled data; 

performing an independent component analysis (ICA) on 
each linearized subset of the sampled data by using the 
data distribution, wherein the ICA generates indepen- 
dent component vectors representing each subset of the 
sampled data; 
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sampling data of known chemical elements representing 
spectral vectors from a database containing the data of 
known chemical elements; 

iteratively performing an act of projecting the independent 
component vectors representative of each subset of the 
sampled data onto the data of the known chemical ele- 
ments, wherein the step of projecting further comprises 
the steps of: 

finding correlation between the independent component 
vectors of each subset of the sampled data and the spec- 
tral vectors of the known chemical elements; 
finding a maximum value of the correlation and a known 
chemical element corresponding to the maximum value; 
creating a list of chemical elements detected on each subset 
of the sampled data, wherein a known chemical element 
corresponding to the maximum value is added to the list 
of chemical elements detected; 
outputting the list of chemical elements detected on each 
subset of the sampled data. 



